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IMTRODUCTICN 


This  final  technical  report  describes  the  research  performed  'jnder  uSAF 
Office  of  Scientific  Research  Contract  F  4462C-74-C-0C23  enritled  "Fracture 
and  Fatigue  of  Bi -Mater i a  I s"  during  the  period  I  December  1973  to  30  Septem¬ 
ber  1978.  The  AFOSR  technical  monitors  were  Mr.  W'lliam  J.  Walker  and 
Lt.  Col.  Joseph  D.  Morgan  III.  \ 

''-S' 

As  suggested  by  the  title  of  J;;his  research  contract,  the  research  is 
directed  toward  a  fundamental  understanding  o''  the  failure  '"echanisms  of  ‘I'he 
filamentary  composites  such  as  boron/aluminum  ana  graoh ita/eoox/  whicn  are 
technologically  the  most  important  examples  of  bi-materials.  At  the  beginning 
of  the  program  there  was  an  attempt  to  use  the  diffusion  bending  crecess  to 
create  a  special  material  which  could  serve  as  a  model  for  bi-material  be¬ 
havior  but  the  diffusion  bonding  process  proved  to  be  much  too  excensive  and 
the  attention  was  thereafter  focussed  on  the  available  filamentarv  corioosite 
mater i a  I s . 

The  penultimate  objectives  of  this  research  are  the  theoretical  founda¬ 
tions  and  the  excerimental  justifica  : orz  for  a  structural  design  merhcdology 
to  be  used  for  the  design  of  filamentary  composite  materials  to  meet  USAF 
damage  tolerance  and  durability  requ  i  rem.ents .  To  that  end  the  work  has  both 
theoretical  and  experimental  phases. 

Ana  I yti ca I  1 y,  the  structural  combination  of  filament  and  matrix  can  be 
viewed  from  rwo  different  geometric  scales.  At  what  is  termed  the  "macrc- 
scooic"  scale,  it  is  assumed  the  filamentary  composite  is  orthotropic  and 
homogeneous.  The  other  scale  is  termed  "micro"  and  views  the  filamentary 
composite  as  a  bi-material. 


On  a  macroscopic  scale,  the  stress  ana  I  vs  is  of  crac:;s  in  ;•  ~'zrrzziz 
materials  and  on  a  microscocic  scale,  -^he  stress  s  i  nga  I  ar  i a"C  stress  dis¬ 
tributions  of  a  crack  lying  at  the  interface  of  a  bi-material  nave  seen 
stud ied. 

The  experimental  work  has  encompassed  the  compression-compression  fatigue 
of  graphite/epoxy  laminates  witn  holes,  1  he  tensile  fracture  of  un i -d i rect Iona  I 
laminates  with  circular  discontinuities,  and  the  tensile  fracture  of  cross- 
plied  laminates  with  various  geometric  configurations  of  discontinuities.  / 

SUMMARY  OF  MAJOR  F IMP  I  MGS 

I.  As  a  result  uf  our  experimental  investigation  of  boron/aluminum  bi- 
materials  we  have  oroposed  a  formula  to  predict  the  rapid  fracture  or 
r+45/O!]^  boron/a  I  umi  num  laminates.  The  motivation  for  this  theory  ccmes 
from  linear  elastic  fracture  mechanics  as  applied  to  homogenous  ma'*‘eria!s 
where  the  fracture  stress,  in  the  presence  of  a  crack  in  a  wide  sheet 

is  given  by 

=  y'TT  (I) 

Thus,  is  a  material  parameter  called  the  fracture  toughness  which  is 

experimentally  determined  and  which  has  dimensions  of  stress  x  length  to 
the  one-half  power.  The  crack  size  effect  (Za  is  the  length  of  the 
crack)  is  embodied  in  the  exoonent  of  "minus  one- ha  If"  which  is  the  order 
of  the  mathematical  singularity  at  the  tip  of  the  crack. 


The  presenT  theory  proposes  an  equation  ter  the  fracture  s'^ress,  c^, 
of  filamentary  composites  of  the  form 

G,  =  H  f2a  )■"’  (2) 

f  CO 

where  H  is  akin  to  but  which  has  dimensions  of  "stress  x  lencth  to 

O  1  ^ 

the  m  power"  (wnich  is  different  than  that  of  The  exoonent  "'"inus 

m"  is  the  order  of  the  singularity  of  a  crack  lying  in  the  matrix  with 
its  tip  at  the  interface  of  matrix  and  filament.  I-;-  has  been  shown  thai' 
the  order  of  the  singularity  is  a  function  of  the  ratio  V|/Uo  cf  the  shear 
moduli  of  the  mamrix  and  filament  and  of  the  two  Possion  ratios  v,  and  v-,. 

Not  much  statistical  data  on  the  shear  moduli  of  boron  filaments,  graphite 
filaments  and  the  epoxy  resins  are  available.  If  ■^he  shear  moduli  of  the 
boron  and  5061  aluminum  are  taken  as  24  million  psi  and  3,3  million  osi 
respectively,  then  the  ratio  of  moduli  for  the  boron/a i umi nun  system  Is 
approximate  I y  .16.  If  the  graphite  filamenx  and  the  ecoxv  resin  tensile 
moduli  are  taken  as  28  and  17  million  psi,  respectively,  then  the  -atio 
is  approximate  I y  .025.  The  most  likely  value  for  m  in  a  boren/a I umi num 
system  is  about  ,33  and  for  a  graph i te/epoxv  system  is  .25. 

The  best  fit  to  our  experimental  data  yields  the  result 

‘^f  -  '0 

^  =  .33(2a  )  ’  (3) 

where  G^  is  the  strength  of  the  laminate  without  holes  or  ether  ^oms  of 


d i sconti nu i t ies . 
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Imbued  with  this  success  we  acplied  this  formula  ro  exoerimerita  I  data 

4 

gatherea  bv  other  experimenters  on  graph i te/epoxy  laminates.  ~h i s  has 
yielded  the  relation 


f  -  23 

=  .43(2a  )  •  ^ 

^0 


for  ro/+4p/90j_  symmetric  lavuDs  and 


a, 

=  .33(2a  ) 

^0 


-.25 


for  C0/+45j_  symmetric  layups, 


(4) 


(5) 


Some  oboervcitions  and  commenrs  on  these  results  may  be  helpful: 

a.  As  predicted  by  the  theory,  the  exponent  in  the  boron/a  I  urn inum 
material  is  larger  than  for  The  graphite/epoxy  material. 

b.  The  experimental  values  are  smaller  than  predicted  by  theory. 

c.  It  should  be  possible  to  correlate  the  factor  H  in  Eg.  (3)  to  the 

stacking  secuence  of  the  layup.  At  any  hole  size,  23,  the  ratio  of 
the  respective  values  for  the  two  graph  i  t-e/epoxy  layups  is 

about  .77  (i.e.,  .33  t  .43).  The  ratio  of  0°  plies  in  the  two 
graph i Te/epoxy  layups  is  .75,  i.e.  one-third  of  the  oLies  in  the 
Q)/;^4^  layuD  is  oriented  at  O"  while  one-fourth  is  the  proportion 
for  the  |^/*45/90^  layuD.  Obviously,  much  more  experimental  worK 

is  recu i red  before  such  a  simple  correlation  can  be  uses. 

d.  The  filaments  in  the  *45  plies  in  a  Layup  contribute  to  the 
"toughness"  whereas  the  ones  at  90"^  do  not  oecause  fracture  can 


occur  solely  in  the  matrix  of  the  90**  plies. 
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e.  Implicit  it  the  crcDOsed  correlaTiv.n  is  the  assL^mction  -hat  cracKs 
in  the  matrix  oerpenc i cu 1 ar  to  the  0°  'i laments  are  the  tritgerinc 
mecnanisn  for  +he  final  failure.  CracKS  in  the  matrix  at  an  angle  to 
the  45°  filaments  have  a  different  kind  of  a  singularity  and  hence 
induce  an  intensity  of  stress  in  these  filaments  insufficient  to 
cause  catastroDhic  failure  in  the  layuD.  However,  the  exoonent  m  'T’av 
be  a  function  of  the  singularities  oresented  by  the  cracks  at  the 
non-zero  degree  filaments. 

2.  Un i -d i rectiona i  t i i amentary/matri x  of  bi-materials  were  tested  in  tne 
dresence  of  circular  hole  and  slit  types  of  discontinuities.^ 

The  failure  process  of  such  a  uni-directional  laminate  begins  as  a  'slit', 
i.e.,  a  crack  propagating  in  the  matrix  oarailel  to  the  filaments.  With 
■^he  aid  of  crack  oroDagation  gages  we  were  able  to  correla'^'e  the  stress 
at  which  the  slit  initiates  against  the  length  of  the  c  i  scont  i  nu  i -f-y . 

This  correlation  leads  us  to  suggest  that  'inear  elasmlc  fracture 
mechanics  may  be  used  and  that  the  ccnset  solit'l'Inc  is  governed  bv  an 
equation  of  the  form 

O  =  H  (23)~'^^  (6) 

s  s 

where  (J_  is  the  far  field  stress  at  the  onset  of  splittinc,  (2a)  is  the 
5 

lencTh  of  the  d i scont i nu i ~v  and  is  a  fracture  touqhness  against 

5 

SD I i tt i ng . 

3.  We  have  tesred  a  small  number  of  ^0/-r45l  and  TOZ-'-cO”!  araoh  i  te/eocxv 
lamina+es  in  the  aresence  of  circular  hole  discontinuities.  There  is  a 


3 


a  difference  'n  the  failure  'nodes  of  these  two  laminates.  Born  sets  c' 
data  do  correlate  with  the  theory  orooosed  in  Section  I  anc  Dreliminarv 
analyses  show  +hat  the  substance  of  commenr  (le)  neeas  additional  inves- 
t i gat  ion . 

4.  Static  compression  tests  of  C'^45/0j^  and  C0/+45l!^  laminates^  were  carried 

out  prior  to  ccmoress ion-ccmoress i on  fatigue  tests,  'hese  laminates, 
which  were  the  facings  of  sandwich  beams,  were  cf  three  di-^^erent  con-fi¬ 
gurations:  (i)  unflawed,  (ii)  a  .25  inch  diameter  hole,  anc  (iii)  a 

.25  inch  diameter  hole  clamoed  by  washers  held  together  by  a  loose  fi'ftinc 
bolt.  The  results  are  shown  in  Fig.  I.  Some  observations  can  be  mace. 

(a)  The  mean  strenqth  of  the  C‘^45/0'’  specimens  is  hicher  than  -for  the 

-  3 

Fo/+45J^  scecimens.  This  agrees  with  the  notion  that  the  -^ree  edge  inter¬ 
laminar  stresses  normal  to  the  olane  of  the  lominare  are  compression  for 
The  iI+45./0j^  system  and  tension  for  the  L0/t45ll^  system.  The  inrer laminar 
tensions  at  the  free  edge  degrade  the  compression  strength. 

(b)  'once  a  noie  is  uni  led,  lue  ■■.eon  otrong+‘'c  -^he  two  svstems  are 

noi  apcreciably  dif-^erent.  This  suggests  that  the  strain  concentrations 
caused  by  the  holes  dominate  the  free  edge  ef-fect. 

(c)  The  clamoed  hole  specimens  exhibit  abc-t  ■^■h.o  co~'  citrengths. 

(d)  The  L'f^5/0j  scecimens  aocear  to  exhibit  a  laruer  scatter  than  do 

-  s 

the  LO/t-lb]  scecimens. 

-  s 

5.  Schema't'ic  recresentations  which  ccmoare  the  sta-f-ic  and  fatigue  failure 
modes  Qf  the  L2;'^5/Gli  and  Il0/+45T^  laminates  as  observed  in  our  -fatigue 
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tests  are  shown  in  Figs.  2 

f o I  I ows : 

(a)  There  are  significanr 
two  laminate  systems. 

(b)  The  nnflawed 

in  both  sTatic  and  fat'coe  '"ai  lores,  ^his  is  '■•ot  *'"'je  'or  ^-e 
system. 

(c)  generally,  the  ■^atigue  '^cdes  o't  ^ailore  are  di''‘ere"~  “re  s^a 

modes  !_w  i +h  the  exceoticns  noted  In  (a)2-  tp  i  5  ocints  ^''e  adce 
com.p  i  ex  i  t  i  es  involved  in  damage  tolerance  and  dorab  i  I  i  a'^al'.-ses. 

In  medals  the  final  failures,  whe~her  cue  "0  o'-'dic  or  s~ati;  adsli- 
catfcns  of  the  loads,  are  the  same  basic  ohencme^’on  and  '-er.ce  can  be 
predicted  pv  the  same  f-ac^jre  mecnanics  formula,  "his  'or'uitcjs 
situation  dees  nc"  acolv  to  the  behavior  o"  ~he  advanced  composites. 
The  varieties  of  ‘a  i  lure  modes  and  the  c  i  ■‘ferences  in  *ai!ure  "'cces 
between  static  and  fatigue  I t^e  acvanced  cerppsi+es  ocirrs  -p  *^e 
need  nor  continued  ressa'"ch. 

(d)  There  dees  appear  *0  be  seme  bene'^iciai  effect  of  ■'■he  clamcinc  acTic^^ 
provided  by  the  washers,  as  can  be  seen  in  .'^ig.  d.  '/ere  da''’a  is  '■e- 
quired  but  if  the  trends  exhibited  in  "io.  d  ore  veri^dd,  then  per¬ 
haps  '"he  beneficial  effect  o^  the  clampinc  provides  pv  mecna^'ical 


esoec^n  i  ve  1  v . 


■erences  r 


specimens  exhibi 


f a  i  i  ^re 


ime  k I  nos 


fas'l'eners  can  be  used  ^or  desicn 


exhi bitea  darace  at  20,000  cvcies.  .he  '”-^^5/0^  crci 

"s  ■' 


of  the  remaining  10  showed  visual  I  v  discernible  "lamace  ar  I',  000  cjc'es. 
The  remaining  soecimen  which  did  nor  exhibit  darrace  a*  10,00  0  c/c'es 

:f  ssec i "ens  exr i - 

5  ■■ 

bited  similar  behavior.  For  the  other  ‘I'wc  crouDS,  !IC/-30l_  anc  _^00/0l_  , 
most  of  the  soecimens  exhibited  damage  at  the  hole  by  3000  cvcies.  ~hese 
-esu  1  ts  lead  '^o  the  "c  I  lowing  com.ments  : 

a.  '-cns  i  derab !  e  camace  Tolerance  was  demcns~r3~ed  bv  tne  ^wc  broucs  .■*".1;'- 
nad  -^‘1?  olies  in  -^haT  ~he  numoers  or  cvcies  *rcm  /isuai  d9~ec~ic' 
damage  am  ■t'he  hole  to  ^inal  ■‘allure  were  large. 

b.  It  can  be  oosrulated  that  the  initial  carnage  at  Tne  hcle  I^me: i a'e ! v 
reduces  the  high  strain  ccncentrati cn  caused  bv  ~he  hole  because  *^.e 
STrain  ccncenTrari  on  exists  cniv  !‘  the  surface  c‘  "‘he  nc  I  e  ana  *'',e 
surrounding  material  memains  conTinuous.  Once  the  S'jrface  c‘  ■‘he  ncie 
is  damagec  bv  spli'*‘s  or  de  I  ami  nat  i  ons ,  then  ■‘he  "soundness''  of  ~‘e 
struct'ure  is  breached  and  the  continuous  strain  ^ieids  whicn  lead  -c 
the  high  concenT^at i cn  cannot  exist,  "Thus,  unlike  '■"era  1  s  and  in 
Oirect  contradistinction  ■t'o  metals.  The  initiaTicn  o'‘  ca^^^ace  a‘‘  a 
hole  leads  to  a  reduction  of  the  strain  CGnceriT'"3‘'i  on  arc  '"3'/  in  carT 


exD  !  a  i  n  the  damage  tolerance  cf  ■^hese  f  i  I  amenta  r’/  systems.  In  metals. 
The  initiation  of  cracks  at  a  hcle  means  a  transition  from  a  stress 
corcen^‘rat i cn  factor  .aDorcach  to  a  fracture  '■"ecnanics  aDcroacn  in 
order  ■‘o  account  ~or  the  theoret  i  ca  I  1  v  in'‘ini'‘e  ccncenf3^  i  a*  ~he 
cracK  tlo.  This  aspect  of  the  behavior  of  comcosite  mamerials,  i 
excloi'‘ed,  would  have  great  technological  imcact. 


c.  These  resu.rs  seem  to  indicate  '^he  ■-anae  of  cvcies  at  -/vnicn  damace 
can  ce  detected  ar  -^he  hole  is  accneciably  le.-s  than  the  ranee  c* 
cycles  wnich  encompasses  railure.  This  is  differen-f  than  for  merals 
because  crack  initiation  in  metals  '■equires  a  scatter  factor  of  four 
but  once  a  crack  has  iiiitiated,  i  fs  growth  to  •'racturo  can  be  encom¬ 
passed  by  a  factor  significantly  le'S  than  four.  Our  excerimenTal 
results  suggest  the  eppesi'^e  for  t  ,e  filamentary  ccmpasi~es.  ~h  i  s 
may  be  due  in  cart  to  *he  ^nef  "^hat  the  observed  damace  mechanisms 
increases  in  comp  I  ex i tv  witn  the  accumulation  of  loaaing  cycles. 

The  ~heoreTical  and  anaivtical  studies  on  the  fracture  and  fatigue  of  bi¬ 
materials  has  addressed  the  macroscopic  stress  analysis  and  the  micro- 
mechanics  of  cracks  (Appendix  A).  These  studies  occur-ed  early  in  -he 
orogram  and  provided  the  theoretical  foudations  for  our  -racture  theory. 
Our  experimental  studies  nave  given  us  a  oh i nomeno I og 1 ca I  Insight  and 
data  base  ucon  which  these  earlier  anaivtical  studies  can  be  aoclied. 

The  stress  distribution  and  smress  intensities  cf  cracks  in  (I)  homoge¬ 
neous  materials  with  complicated  geemetry  (2)  bi-material  in+erface 
(3)  anisotropic  materials  have  been  made  amenable  top  simoiified  ana 
economical  solution  by  the  creation  of  crack-containing  finite  elements. 

A  familv  of  bi-material,  homogeneous,  and  anistreoic  soecial  crack 
elemen+5  in  plane  and  axi symmetric  cases  have  been  develoced.  These 


i  nc  1  ude 
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3.  Hcmoceneous  mixed-'^cde  olane  cr^ack  ele'^enr 

b .  Bl-material  debonding  olane  crack  ele'^ent 

c.  Bi-'raterial  peroend  i  cu  I  ar  plane  crack  elerenr 

d.  Axi symmetric  penny-shaped  crack  element 

e.  Ax [symmetric  circumferential  crack  elemenr 

*.  AxisvmmeTric  bi-marerial  debcnding  crack  eie^enr 
c.  Anisotrocic  olane  crack  elemenr 
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APPEMOrX  A 

FRACTURE  MECHANICS  OF  COMPOSITES 

Introduction 

The  concepts  of  fracture  mechanics  are  now  being, applied  to 
the  study  of  crack  extension  phenomena,  both  slow  and  catastrophic, 
in  isotropic  homogeneous  materials,  and  is  of  great  importance 
wherever  high  strength  metallic  alloys  are  used.  Based  on  energy 
equilibrium  considerations,  Griffith  (1920)  first  developed  a 
fracture  theory  for  brittle  materials  such  as  glass  [4],  Griffith's 
theory  postulates  that  fracture  occurs  when  the  decrease  of  elas¬ 
tic  strain  energy  per  increment  of  crack  area  is  equal  to  the 
increase  of  surface  energy  over  the  same  area.  However,  the 
surface  energy  -^or  a  solid  is  difficult  to  calculate  and  measure. 

To  overcome  this  difficulty,  Irwin  introduced  the  concept  of 
"stress  intensity  factor"  K,  which  is  related  to  the  amplitude 
of  the  stress  field  near  the  tip  of  a  crack  [5] .  In  fact,  it  can 

be  shown  that  the  stress  intensity  factor  is  related  to  the  rate 

3  U 

of  energy  release  per  unit  area,  G  =  and  hence  to  Griffith's 

surface  energy.  The  Irwin-Griffith  fracture  criterion  states 

that  catastrophic  crack  propagation  occurs  when  G  or  K  reaches  a 

critical  value  G  or  K  and  this  is  regarded  as  a  material  prop- 

c  o 

erty,  called  "the  fracture  toughness"  of  the  material. 

It  can  be  shown  that  under  general  loadings,  the  stress  dis¬ 
tribution  of  an  elastic  body  in  the  presence  of  cracks  can  be 
superimposed  from  three  separate  cases:  "mode  1  which  is  the 
opening  mode",  "mode  2  which  is  the  shearing  mode",  and  "mode  3 
which  is  the  out-of -plane  shearing  mode"  [6] .  For  homogeneous 
isotropic  materials,  values  and  formulae  of  stress  intensity 
factors  obtained  from  various  methods  of  analyses  can  be  found 
in  various  handbooks  [7].  A  large  body  of  laboratory  data  has 
also  been  generated  and  there  is  sufficient  evidence  that  frac¬ 
ture  toughness  is  a  material  constant,  independent  of  crack 

c 

length  and  crack  configuration. 
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Figure  7.  Accuracy  Test  of  Present  Finite  Element  Solutions 


Figure  8.  Stress  Intensity  Factors  of  a  Center  Crack  with 
Boron/Aluminum  Laminates  Properties 
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Figure  9.  Errors  in  Approximating  Isotropic  Stress  Intensity 
Factors  for  a  Center  Crack  with  Unidirectional 
Graphite/Epoxy  Properties 
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Figure  10.  Stress  Intensity  Factors  of  a  Center  Crack  with 
Graphite/Epoxy  Properties 
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Figure  11.  Stress  Intensity  Factors  of  a  Double-Edged 
Crack  with  Boron/Aluminum  Properties 
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re  13.  Errors  in  Approximating  the  Isotronic  Stress 
Intensity  Factors  for  a  Double-Edge  Crack 
With  Unidirectional  Graphite/Epoxy  Properties 
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A  Debonding  Crack  Along  the  Bi-material  Interface 


constant 


ooJO/ 


J  / 


(a)  Geometry 


(b)  Mesh  Plan 


Figure  17.  A  Crack  Near  a  Bi-material  Interface 
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Figure  18.  Comparison  of  Finite  Element  Solutions  with 
Other  Solutions 
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(a)  Geometry 


(b)  Mesh  Plan 


Figure  19.  A  Crack  in  a  Softener  Strip 
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(b)  Mesh  Plane 


Model  for  Axisyrrunetric 


re  22.  Finite  Element  Solutions  of  a  Penny-Shaped  Crack 


re  23.  Finite  Element  Solutions  of  an  Anisymmetric 
Circumferential  Crack 


ble  1,  Roots  X  of  Characteristic  equation  (8) — to  be  Used  for  Finite  Element  Analysis 
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_Table  2  Accuracy  of  Finite  Element  Solutic 


1 

1 

a/B 

H* 

f.e.m. 

Bowie  Solution 
[477  3 

.10 

1.00 

1.01 

.20 

1.04 

1.05 

.30 

1.10 

1.11 

.  40 

1.18 

1.19 

.50 

1.28 

1.29 

.60 

1.  40 

1.41 

.70 

1.58 

1.58 

.  80 

1.86 

1.  85 

*  1.  H  =  Kj/c/a 

2.  =  1,  B2  =  .32,  L/B  =  2.0 


Table  3.  Stress  Intensity  Factors  and  for  a  Crack  Along 
the  Bi-Material  Interface 

=  1  psi/  Vj^  =  0.30/  ~  ^  psi/  ^  psi 

20  X  20-inch  panel,  crack  length  =  2  in.,  plane  stress 


V1/V2 

^2 

(psi) 

^1 

^2 

EXACT* 

F.E.M. ** 

EXACT* 

F.E.M. 

1 

.3 

1.00 

1.000 

■  1.01 

0.0 

0.0 

3 

.3 

.53 

.988 

.997 

.0724 

.0728 

10 

.3 

.37 

.968 

.975 

.1171 

.1176 

23.1 

.35 

.38 

.966 

.972 

.1208 

.1213 

100 

.3 

.31 

.953 

.959 

.1391 

.1395 

144.2 

.35 

.36 

.960 

.965 

.1303 

.  1306 

1,000 

.3 

.30 

.952 

.957 

.1415 

.1418 

♦Solution  for  an  infinite  panel, from  Rice  and  Sih  125] 
♦♦Finite  element  solution. 
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Table  4.  ENERGY  RELEASE  RATE  G  FOR  A  CRACK  ALONG  THE 
BI-MATERIAL  INTERFACE 


=  22.5,  =  0.2,  Vj  =  0.35,  0^^  =  1  psi,  =  1  psi 

^xx2  ~  psi/  20  X  20-inch  panel,  plane  stress 


Half 
Crack 
Length 
(in. ) 

Half-panel 

b  ujL  din 

Energy 

G  = 

dU/dc 

k 

2 

EXACT* 

F.E.M.** 

EXACT* 

F.E.M. 

EXACT* 

F.E.M. 

0.80 

969.1628 

25.04 

— 

0.862 

0.867 

0.125 

0.122 

0.81 

969.4165 

25.35 

25.37 

0.868 

0.872 

0.124 

0.122 

0.82 

969.6734 

25.67 

25.69 

0.873 

0.878 

0.124 

0.122 

0.84 

970.1970 

26.29 

26.18 

0.884 

0.889 

0.124 

0.122 

0.86 

970.7337 

26.92 

26.84 

0.895 

0.900 

0.123 

0.121 

♦From  Reference  1 4  8] 

**AU/AC  calculated  by  taking  differences  between  half  crack  length 

and  half-panel  strain  energy,  vising  consecutive  lines  in  ^he  table. 


Table  5.  ENERGY  RELEASE  RATE  G  FOR  HOMOGENEOUS  MATERIAL 

E  =  1  psi,  =  1  psi,  v  =  0.30,  20  x  20-inch  panel, 

plane  stress 


Half 
Crack 
Length 
(in. ) 

Half-panel 

Strain 

Energy 

G  = 

dU/dc 

k 

1 

EXACT 

F.E.M. 

EXACT 

F.E.M. 

0.80 

141.0115 

2.5T. 

— 

.894 

.900 

0.81 

141.0370 

2.55 

2.55 

.900 

.906 

0.82 

141.0629 

2.58 

2.59 

.906 

.911 

0.84 

141.1156 

2.64 

2.64 

.917 

.923 

0.86 

141.1697 

2.70 

2.71 

.927 

.934 

Table  6.  Mode  1  Stress  Intensity  Factors  fe  and  fe  for  a  Crack  Normal  to  the 
Bi-Material  Interface  °  " 
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APPENDIX  A, 


Stress  distribution  near  the  bi-material  crack  tip 
I.  Crack  lying  along  the  interface 

For  +  iXj,  n  =  1,  2,  3  , - 

^=Z1?e[A„t  b„fo,+2}te  — X^-"  -(.K-Oie,  "  ] 

Tl 

«^y  =Z1?e  {A„r CA'Oi 


_  /O  f  ^  f'  \  -aCA-OS  ZtArU0 

<Zy=Z  1?e  l  Afl'T  (xbJ[-\•^e  -/-e-  -{-CX-ofe 


For  X^=n,  n=  1,  2#  3^ 


=ZA„r^'^|  CCA+3)C^^(A-Oe-(VOc^56\-3;0jCn 

4  [  C  A+ 1 J  CA„- 1 )  0  -  U5i«-  (A„'3;  ej  j 

“Z7),ir  ^  I  lrCA„- I  )cosC^-\)9-^ 

r-(A„'2;  s-<^0^.-ue  N-CA J  1 

5  {  C-'^Z+D  (7'„-'U5^0'3J)^(;„ 

+  ( (>„-(;  ‘^^c^-l)0-c;^,-lJc^sc\-3JsJd„\ 


where 


1-1 , 1-.^  , 

'  *•!  .  1 


f 


in  material  1 


in  material  2 
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II.  Crack  normal  to  the  interface 


<3x  =  Z  ]_>„  r"'  b„  [ 

-  CA.-O  L-{'g,coiLK-3)e-j^s^  j 

+  (A„'U('fRw(>„'3;0-'ft5XA„'3j)5}j'' 

{zr'"  t5,  r  + 


where 


in  material  1 


(3 

C> n  ) -  -  2 ^  _  £7<  K^]  /v  (\  ) 

JPC^ri)^  2c< (c<'f'c<^J)<:^5✓^ffT^— 


and  in  material  2 
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APPENDIX  A. 2 

Stresses  and  Displacements  Used  for  the 
Derivation  of  Axisynmetric  Crack  Elements 


%  ^2  \  f  >  [(3  'A-  iO  0^-2;  9  -f  cv  z;  CA- 


<!-'  1.  V 


t. 
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A^  =  y  f  v  . 
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where 

(1)  X  =  I  /  I  /  J  '  •  •  •  •  '  =  2,  Y  =  -1 

C  "f"  rrt  CA-  3j)  "A. 

^  CA  "23  "XH  6A"i3  —A 

-^4^  2f  yn -H  5!a''4"' C)^-+A3 
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-^1  =  (A'2-jV  I^'Q  "A 
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APPENDIX  A. 3 

Relation  Between  the  Stress  Intensity  Factors 
^1'  ^2  Constants  bj 


j  j  j+n  ,  b j ,  real  constants.  n  =  No. 

stress  (or  displacement)  terms  assumed. 


1.  Homogeneous  Isotropic  Materials 

=vtK 

2.  Homogeneous  Anisotropic  Materials 

=  yi:  (Sz-s,V-5i  'o' 

3.  A  Debond  Crack  Along  the  Bi-Material  Interface 

4.  A  Crack  Normal  to  the  Bi-Material  Interface 

is  defined  in  Appendix  2. A. 


Penny-Shaped  or  Circumferential  Crack 


'  /r-  ” 


of 


\J^ 

APPENDIX  A. 4 

Material  Properties  for  the  Fini 
Element  Analysis 


I.  Boron/Aluminum 
(i)  (0)  Layup: 

E  =  20.1  msi,  E 

X  2 

G  =7  msi,  S,  = 

xy  1 

(ii)  (±45/02)  Layup: 

E  =  20.3  msi,  E 
X  y 

G  =  8.7  msi,  S, 

xy  X 

(iii)  (90)  Layup: 

E  =  33.1  msi,  E 
X  y 

G  =7  msi,  S,  = 

xy  1 


=  33.1  msi, 
1.5180i,  S2  =  . 

=  26.5  msi, 

=  l-0848i,  S2  = 

=  20.1  msi, 
1.9480i,  S2  ==  . 


II.  GranH-’ te/Epoxy 
(i)  (0)  Layup: 

E^  =  1.7  msi,  Ey  =  17  msi,  = 

G  =  0.65  msi,  S,  =  1.5918i,  S- 
xy  '1  2 

(ii)  (90)  Layup: 

E^  =  17  msi,  E„  =  1.7  msi,  v  = 

X  y  y  A 

G  =  .65  msi,  S,  =  5.0338i,  S_  = 
xy  1  2 


=  .25 
5134i 

=  .33 
.  8068i 

=  .15 
6588i 


.21 

=  .198661 

.  021 


6282i 
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Appendix  f 
COUPLING  ACTIVITIES 

The  Principal  Investigator,  Prof.  James  W.  Mar,  was  very  active  during 
the  years  of  this  contract  In  the  USAF  Scientific  Advisory  Board  (SAB).  As  a 
member  of  the  SAB  he  was  assigned  at  various  times  to  the  Division  Advisory 
Grouos  of  the  Aeronaurical  Systems  Division,  the  Armament  Develoomen-r  ana 
Test  Center,  the  Foreign  Technology  Division,  the  Logistics  Cross  Panel,  and 
the  Science  and  Technology  Advisory  Group.  Some  of  the  advisory  activities 
which  were  impacted  by  the  research  activixies  of  this  contract  are  as 
fo I  I cws ; 

1.  Chairman  of  the  USAF  SAB  Ad  Hoc  Committee  on 

Advanced  Comoos ites  Technology,  1976. 

2.  Chairman  of  the  USAF  SAB  Ad  Hoc  Committee  on 

the  B- I  Strucutre,  1973-1977. 

3.  Che '"man  of  the  USAF  AFSC  STAG  Ad  Hoc  Com.mittee 

on  3-0  Carbon/Carbon  Nozzles,  1977. 

The  Principal  Investigator  during  this  period  was  also  a  memcer  cf  "he 
following  advisory  Committees: 

1.  NASA  Research  and  Technology  Advisory  Committee 

on  Materials  and  Structures, 

2.  NASA  Space  Systems  Committee. 

3.  OOD  Defense  Science  Board  Task  Force  on  V/STCL. 

4.  AGARD  Structures  and  Materials  Panel. 

Of  special  note,  three  students  who  participated  in  the  research 
activities  of  this  contract  are  now  working  in  the  field  of  advanced 
composites.  Their  present  assignments  are  due  in  large  measure  to  their 
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involvement  with  their  research.  They  are  listed  as  follows 
Dr.  K.  Y.  Lin,  Boeing  Commercial  Airolane  Co. 
i^r.  David  Maass,  Sikorsky  Aircraft 
2nd  Lt.  Karyn  Knoll,  AFML 


